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ABSTRACT 

We present a radiative transfer model for the circumstellar dust shell around a Large Magellanic Cloud (LMC) long-period variable 
(LPV) previously studied as part of the Optical Gravitational Lensing Experiment (OGLE) survey of the LMC. OGLE LMC LPV 
28579 (SAGE J05 1306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have Spitzer broadband 
photometry and spectra from the SAGE and SAGE-Spec programs along with broadband UBVIJHK S photometry. By modeling this 
source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce 
the spectral energy distribution of the source using a mixture of amorphous carbon and silicon carbide with 15% SiC by mass. The 
grain sizes are distributed according to the KMH model, with y = 3.5, a mm = 0.01 jim and «o = 1-0 jim. The best-fit model produces 
an optical depth of 0.28 for the dust shell at the peak of the SiC feature (1 1.3 yum), with an inner radius of about 1430 R Q or 4.4 times 
the stellar radius. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s -1 , we obtain a dust 
mass-loss rate of 2.5 x 10~ 9 M Q yr~' . We calculate a 15% variation in this mass-loss rate by testing the sensitivity of the fit to variation 
in the input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 
13.7 yum feature seen in the spectrum. We find that a combination of CO and C 2 H 2 gas at an excitation temperature of about 1000 K 
and column densities of 3 x 10 21 cirT 2 and 10' 9 cirT 2 respectively are able to reproduce the observations. Given that the excitation 
temperature is close to the temperature of the dust at the inner radius, most of the molecular contribution probably arises from this 
region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L Q . For an effective temperature of about 3000 
K, this implies a stellar mass of 1.5-2 M Q and an age of 1-2.5 Gyr for OGLE LMC LPV 28579. We calculate a gas mass-loss rate of 
5.0 x 10~ 7 M yr~' assuming a gas:dust ratio of 200. This number is comparable to the gas mass-loss rates estimated from the period, 
color and 8 fim flux of the source. 

Key words, circumstellar matter, infrared: stars, stars: asymptotic giant branch 
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?3 i 1. Introduction The large numbers of AGB stars makes them significant con- 

, , . __, _ , . tributors to the chemical evolution of galaxies. In particular, 

The asymptotic giant branch (AGB) is the final stage in the c . rich AGfi starf , arg ±e a ^ major sourcef , of carbon 

evolution of low- and intermediate-mass stars (0.8 to 8 M ). atQms (Dwek 199g) afld carbonaceous dust ins Jp-^li998: 

Products of the nuclear reactions in the AGB star interiors are Matsuura et al 2009) in a al It is therefore cmcial tQ deter . 

brought to the outer regions where they form molecules and, far- mine ^ osition of AGB star dust and the total injection 

ther from the star, dust grains. Carbon-rich (C-r ich) AGB stars rate of Ms dust ^ ^ ISM 

are cr eated by the process of the third dredge-up (Iben & Renzini 

1983). Slow winds, with typical expansion velocities of 10 km The stud y of Galactic carbon stars is complicated due 
s" 1 , feed the processed material (gas and dust) back into the sur- to the hl 8 h ( and often unknown) line-of-sight extinction to- 
rounding interstellar medium (ISM) at rates up to 10" 5 - 10" 4 wards these sources - Des P lte thls fact ' a lot of information has 
M yr-\ where it is eventually locked into newly forming stars. been 8 ained about the circumstellar dust around these stars. 

The Large Magellanic Cloud offers an exceptional environment 

* Figures 3 and 4 ar e only available in electronic form via for the study of the mass-loss return and enrichment process 

http://www.edpsciences.org by AGB stars. Recent near-infrared (NIR) and mid-infrared 
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(MIR) surveys (e.g. , 2MASS, ISkrutskie et all d2006l): D ENIS. 
lEpchtein et all (11994 and SAGE. iMeixner et al.l (l2006h ) have 
revealed tens of thousands of AGB star candidates in L arge 
Magellanic Cloud (LMC). The SAGE dMeixner et alJl2006h sur- 
vey has observed the LM C in the Spitzer Space Telescope 
(Spitzer; iGehrz et al.ll2007h Infrared Array Camera (IRAC; 3.6, 
4.5, 5.8 and 8.0 pm) and Multiband Imagi ng Photometer for 
Spitz er (MIPS; 24, 70 and 160 pm) bands. iBlum et all d2006l) 
and I Srinivasan et aTl d2009l) (hereafter, Paper I) used SAGE pho- 
tometry to identify about 7, 000 carbon-rich AGB star candi- 
dates in the LMC. Follow-up spectroscopy using the Infrared 
Spectrograph (IRS) on Spitzer has also been per formed on se- 
lected SAGE targets as part of the S AGE-Spec ( Kemper et al. 
program. 

Detailed radiative transfer (RT) studies conducted on 
Galactic and extragalactic carbon stars in the past have helped 
identify many dust species as possible constituents of their cir- 
cumstellar envelopes. Here we concentrate on two species that 
significantly affect the near- and mid-IR photometry and spec- 
troscopy, carbon dust (amorphous car bon or graphite) and sili- 
con c arbide (a or ft type). IRC+10 216 (Neu gebauer & Leig hton 
1969) is the best-studied carbon star in the galaxy, and a 
majority of early st udies involved the modeling of its cir - 
cumstellar shell (e.g.. Hones & Merrilllll976l; iFazio et alj [l980). 
iRowan-Robinson & Harr is (1983) used their RT code to fit 
the circumstellar envelopes (CSEs) of 44 Galactic carbon stars 
and found that amorphous carbon dust (hereafter AmC) pro- 
du ced better fits than graphite. A similar result was ob tained 
bv lMartin & Rogers! dl987l) . ISkinner & Whitmorel dl988l) found 
that the strength of the SiC feature in carbon stars was in- 
versely proportional to the total mass-loss rate which they 
explained as being due to a decrease in the SiC/AmC ra- 
tio with increasing mass-loss rate, thus decreasing the promi- 
nence of the SiC feature in comparison to the AmC contin- 
uum. This decrease in SiC/AmC ra t io wi th increasing mass 
loss was verified by IChan & Kwokl d 19901) and explained as 
being due t o carbon star evolution. Similar trends w ere also 
reported by Lorenz-Martins & Lefevre ( 19931 11994 and by 
IGroenewegenl d!995l). IGroenewegenl d 19951) . ISpeck et all dl997l) 
and iBlanco et alJ dl998l) compared the results of models using 
a- and /?-SiC, finding that the dust containing g -SiC was able 
to rep roduce the spectra of most of their sources. iBressan et alJ 

(1998) included a treatment of AGB dust shells in their stellar 
population models, and they generated isochron es for compari- 
son w ith IRAS data for Miras and OH/IR stars. Ivan Loon et alJ 

(1999) performed RT modeling on the ISO spectroscopy of LMC 
AGB stars to identify the chemistry of the circu mstellar du st and 
calculated their mass-loss rates and luminosities. Suh (2000) ob- 
tained empirical opacity functions for amorphous carbon dust 
based on their RT mo deling of IR spectra a s well as laboratory- 
measured optical data. Groenewegen ( 2006) presented synthetic 
photometry for O-rich and C-rich AGB stars from the results 
of RT models that spanned the relevant range of stellar and 
dust shell parameters seen in Galactic AGB stars. In particular, 
they considered two kinds of dust species for carbon star dust 
shells: AmC and a mix ture of AmC and 15 % a -SiC by mass. 
The theoretical work of Mattsson et alJ (2008]) and lWachter eTal] 
(2008) showed that the mass-loss rates of carbon stars may not 
be sensitive to metallicity. Numerous recent studies using Spitzer 
spectra o f AGB stars in low-metallicity Loca l Group galax- 
ies (e.g.. ISloan et alJ )200d jZiilstra et al. 2006; lLagadec et alJ 
120071: iMatsuura et alJl2007t ISloan et al.ll2008h seem to support 
this cl aim. Similar results were obtained bv lGroe newegen et alJ 
(2007), who found in their modeling of the IRS spectra of 60 



carbon stars in the Magellanic Clouds that the trend of mass- 
loss rates with p eriod of luminosity was co mparable to that of 
Galactic sources. iGroenewegen et al] (12009) extended this study 
to a larger sample, including 110 carbon stars, and compared 
their mass-loss rates and luminosities with those predicted by 
evolutionary models as well as dynamical wind models. 

Our aim is to develop a grid of RT models for dust shells 
around carbon stars so that we may fit the SAGE photome- 
try of carbon star candidates and derive their mass-loss rates. 
These mass-loss rates will enable the calculation of the total 
carbon-star mass-loss return to the LMC. As a first step, in this 
paper, we determine a representative set of dust properties for 
LMC C-rich AGB stars to be used as input to the modeling. 
We present a 2Dust radiative transfer model for the circumstel- 
lar shell around the variable carbon star OGLE LMC LPV 28579 
(SAGE J051306.40-690946.3) with the primary goal of deriv- 
ing the dust properties with a physically realistic model of the 
source, as constrained by available data. We have also used this 
approach in Sargent et al. 2010 (in press) for the O-rich AGB 
stars and red supergiants. This paper is organized as follows. 
We describe the observational data for OGLE LMC LPV 28579 
(hereafter LPV 28579) from various studies in Sect. [2] In Sect. [3] 
we provide details of our radiative transfer model for the circum- 
stellar dust, as well as a simple model for the molecular features 
observed in the spectrum. We discuss the results of the models 
in Sect. [4] 



2. Observations 

The choice of LPV 28579 for this study was motivated by the 
availability of both SAGE photometry and S AGE-Spec spectro- 
scopic data. The SAGE-Spec program selected a bright subsam- 
ple of the SAGE AGB candidates for good quality spectra. This 
requirement means that the spectroscopic sample is biased to- 
wards redder colors, and that LPV 28579 is redder than most 
carbon stars in the LMC. Nevertheless, it is only moderately op- 
tically thick, and exhibits features found in typical carbon stars. 
Its spectrum (see Sect. I2.3l for details) shows the overall contin- 
uum emission from carbon dust as well as a significant 1 1.3 pm 
feature due to SiC, which makes it a good testbed for the dust 
properties around LMC carbon stars. IGroenewegenl d2004) clas- 
sified the star as an obscured AGB candidate, whil e in Paper I 
we lab elled it an "extreme" AGB star based on the IBlum et aTl 
(2006) selection criterion ( J - [3.6] > 3.1 mag). Examination of 
the SAGE-Spec spectrum (Section l2.3t shows the 11.3 pm sili- 
con carbide feature to be in emission. Moreover, results from the 
present modeling work (Section [3.1.5l ) suggest at best a moder- 
ate 1 1 .3 pm optical depth. Based on this information, LPV 28579 
would not be considered an extreme carbon star dSpeck et al.l 
120091) . 

Photometry for LPV 28579 is available from many recent 



LMC surv eys (e . g.. IZebrun et al l 1200 It lEpchtein et alJ 11 1999: 
ICutri et al] 120031: iMeixner et all 120061: iKato et al] 120071) . en- 
abling us to constrain its spectral energy distribution (SED). We 
have combined the SAGE data with photometry fr om the optical 
Mage llanic Clouds Photometric Survey (MCPS; Zaritsky et al. 
1 19971) as well as t he 2 micron All Sky Survey (2MASS, 
Skru tskie et al.ll2006l) . We discuss the optical and NIR variabil- 
ity observations in Sect. 12.11 the SAGE photometry in Sect. 12.21 
and the SAGE-Spec data in Sect. 12.31 In Fig. Q~l we show the 
full range of photometry and spectroscopy values measured to 
illustrate the source variability. In this paper, we only model the 
SAGE Epoch 1 photometry, using the optical and NIR data along 
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with the IRS spectrum to constrain the shape of the resulting 
SED. 

2.1. Variability in Optical and NIR Data 

LPV 28579 was observed as par t of the Optical G ravitational 
Lensing Experiment (OGLE-II; Zebru n et al.l I2QQ lT) survey of 
the Magellanic Clouds. The / band light curve is available as 
part of the OGLE-II variable star catalog (Szvmahski 2005b as 
well as th e OGLE-III list of long -period variables (LPVs) as de- 
scribed in lSoszvnski et alj d2009l) . Based on its light curve , LPV 
28579 was classified as a Mira-type LPV. Ilta et al] (I2004 and 
iGroenewegenl d2004l) crossmatched th e OGLE-I I data with the 
IRSF LMC Survey dKato et all 1 20071). DENIS dEpchtein et all 
1999) and 2MASS All-Sky Release dCutri et al.l2003l) NIR cata- 
logs and fit the light cur ves t o obtain variability p eriods and am- 
plitudes. llta et al.l(l2004l) and lGroenewegenl d2004l) calculated the 
period of varia tion to be 361 d and 358.6±0.136 d respectively. 
Ilta et alJ (|2004) also reported a mean / magnitude of 17.483 mag 
and a peak-to-peak amplitude of 1.865 mag. The OGLE-III cat- 
alog lists the mean magnitude as 17.586 mag. The light curve 
available from the OGLE-III websitefl shows that LPV 28579 is 
multi-periodic, with a primary period of 356.2 d (peak-to-peak 
amplitude 1.587 mag) modulated by a slow variation of 5747 d 
(peak-to-peak amplitude 0.784 mag). Like many of our extreme 
A GB candidates identi fied in Paper I, LPV 28579 is absent from 
the lFraser etaT1d2008l) MACHO catalog of LMC LPVs. This is 
most likely due to the moderately high circumstellar extinction 
at optical wavelengths. The primary period for LPV 28579 is 
very close to those of the sources in the "one year ar tifact" list, 
but the culling of these artifacts in lFraser et al.l d2008l) was done 
in a manner as to avoid accidental removal of genuine sources. 
Based on its brightness (7<f s =10.9 mag) and period, LPV 28579 
wo uld probably fall o n period sequence 1, the sequence found 
bv lWood et al.l d 19991) to be predominantly populated by Miras 
in their fundamental mode of pulsation. While we do not have 
JHK S light curves, the 2MASS and IRSF photometry for LPV 
28579 capture some of the NIR flux variatiorQ (Fig. [I}. 

The optical properties of LPV 28579 are fairly normal for 
LMC carbon stars - the pulsation period of ~ 1 yr is average for 
typical LMC carbon Miras. It has a relatively small amplitude in 
comparison to other carbon Miras in the OGLE-III catalog (am- 
plitudes ranging from 0.8 mag to 5.7 mag), and this points to a 
somewhat hotter central star than a very cool one. The moderate 
amplitude combined with the ~ 1 yr period make OGLE LMC 
LPV 28579 a typical member of the LMC carbon Mira family. 
Furthermore, while LPV 28579 has very red colors (J-K-3.27 
mag) , it is nowhere nearly a s red as, e.g., IRC+10216 (J-K=6.15 
mag: IWhitelock et al.ll2006h . 

2.2. SAGE Photometry 

The SAGE survey dMeixner et al.ll2006l) imaged a ~ 7° x 7° area 
of the Large Magella nic Cloud (LMC) with the Spitzer Infrared 
Array Camera (IRAC: lFazio et al.l2004l 3.6, 4 .5, 5.8 and 2, urn) 
and Multiband Imaging Photometer (MIPS: iRieke et al.| [2004. 
24, 70 and 160 fim). Two epochs of observations separated by 



1 http://ogledb. astrouw. edu.pl/~ogle/CVS/ 

2 The difference in A" s magnitu de between the two ep ochs is 0.5 mag. 
Comparing this result to Fig. 3 of Whitelock et al. ( 2003), which shows 
the A^-band pulsation amplitude as a function of period, we find that 
the amplitude corresponding to the primary period for LPV 28579 is 
-0.4-0.8 mag. 



three months were obtained to constrain source variability. IRAC 
and MIPS epoch 1 and epoch 2 archive photometry for LPV 
28579 is available as part of SAGE data delivered to the Spitzer 
Science C enteiEl (SSC). By comparing data from both epochs 
of SAGE, Viih et all d2009l) identified infrared variable sources 
in the LMC of which -81% are AGB candidates. L PV 28579 
is one of the ~2 000 SAGE variables discovered by IViih et all 
(2009). FigureQ]shows both epochs of SAGE data. In this paper, 
we model the Epoch 1 data (the brighter set of fluxes in Fig.[T|) 
and use the Epoch 2 photometry only to estimate the range of 
variability of the source (see Sect. 13. 1 .Tb . 

2.3. SAGE-Spec Data 

Spectroscopic follow-up observations to SAGE using the 
Infrared Spectrograph (IRS) aboard Spitz er were performed a s 
part of the SAGE-Spectroscopy survey (Kem per et alJ 12010). 
The spectra were reduced using tech niques described i n the 
SAGE-Spec Data Delivery Handbook dWoods et alJl2009l) and 
summarized briefly as follows. LPV 28579 was observed in both 
spectral orders of the Short-Low (SL; 5.2-14.3 pm, R-60-120) 
and Long-Low (LL; 14.3-13.7 /mi, R~60-120) instruments on 
board Spitzer at two nod positions. Data from the S15.3 and 
S17.2 pipelines for SL and LL, respectively, were obtained from 
the SSC for LPV 28579 (AOR # 22415360). SL sky subtrac- 
tion was performed by subtracting the observations in one order 
from the other while keeping the nod position fixed, while LL 
sky subtraction involved subtracting one nod position from the 
other while keeping the order fixed. Bad pixels on the detector 
arrays were replaced with values determined by comparing to 
the local point-spread function^. 

The width of the apertures used for point source signal ex- 
traction increased in proportion to the wavelength; this involved 
using the profile, ridge and extract module^- Flux densities were 
obtained by calibrating the raw signal extractions using observa- 
tions of the standard stars HR 6348 (K0 III), HD 166780 (K4 
III) and HD 17351 1 (K5 III). While all three standard stars were 
used for LL, only the first was used for SL. Where available, 
the mean flux from multiple data collection events (DCEs) for 
a single module/order/nod was computed, with the standard de- 
viation from the mean providing a measure of the correspond- 
ing uncertainty. Spikes not eliminated by replacement of bad 
pixel values were effectively removed by replacing any uncer- 
tainty larger than the average in the local wavelength range by 
the uncertainty from the nod position closer to the local average. 
Bonus-order spectra (corresponding to a small segment of the 
first-order spectrum obtained by the second-order slit) were av- 
eraged with spectra from the second and first orders over wave- 
lengths they shared in common. 

The IRS spectrum for LPV 28579 (Fig. HJ shows a broad 
1 1.3 pm SiC emission feature combined with a featureless con- 
tinuum characteristic of carbonaceous dust. LPV 28579 does 
not show any particular signs of the 30 micron fe ature which is 
strong in the more extreme Galactic carbon stars dHrivnak et al.l 
2009) and is also seen in very red LMC carbon stars in the 



3 Data version S13 available on the SSC website, 
http : //ssc . spitzer . caltech . edu/legacy/sagehistory . html 

4 Handbook version 1 (30 July 2009) available on the SSC website, 
http : / /data . spitzer . caltech . edu/popular/sage-spec/ 

5 For details, please refer to the IRS Data Handbook: 
http : //ssc . spitzer . caltech . edu/IRS/dh/ dh32 . pdf 

6 These modules are available in SPICE, the Spitzer IRS Custom 
Extraction package. 
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SAGE-Spec sample. There is also a 13.7 /im C2H2 feature. The 
strength of this feature cannot be reproduced by a model photo- 
sphere alone, pointing to its origin in the extended atmosphere 
and circumstellar shell (see Matsuura et al.l l2006h . We model 
both the circumstellar dust as well as the molecular gas in the 
extended atmosphere in Sect. [3] We note two trends at longer 
wavelengths: there is a step-down in the flux at ~20.5 //m, and 
the flux errors for the entire Long-Low 1st order (LL1) wave- 
length range are on average very large. These are suggestive of 
sky-subtraction issues, so we checked the raw data for LL1. We 
find that there is a noticeable amount of extended emission in 
LL1 due to an extremely red source creeping into the slit, which 
affects the background subtraction. 

The optical, NIR and MIR photometry as well as MIR spec- 
troscopy described above are incorporated into Fig.Q] The two- 
epoch photometry shows that there is no significant change in 
the shape of the MIR SED. We have taken advantage of this fact 
and scaled the SAGE-Spec spectrum down to the SAGE Epoch 1 
5.8 fj.m flux to enable easy comparison to photometry. This scal- 
ing changes the flux predicted for the spectrum in the 5.8 fim 
band by -5%. 

3. Analysis 

3.1. 2Dust Radiative Transfer Models 

We model the emission from LPV 28579 using the two- 
dimen sional radiative transfer code 2Dust dUeta & Meixnerl 
I2003I) . For simplicity, we assume a spherical geometry for this 
study. 2Dust is capable of radiative transfer of axisymmetric sys- 
tems. This may be needed for describing the most evolved AGB 
stars as well as their post-AGB counterparts in the model grid. 
Although AGB mass loss appears to be largely spherical, the 
successors to these objects, the post-AGB stars, are expected to 
be bipolar (as are their Galactic counterparts). The dust around 
post-AGB stars is similar to those around the most evolved AGB 
stars, therefore we allow for the model grid to cover post-AGB 
stars and pre-planetary nebulae (PPNe). 

The goal of the present study is to extract a set of dust proper- 
ties for LPV 28579 that will also provide a range of reasonable 
properties for our entire C-rich AGB sample. As input, 2Dust 
needs information about the central star (effective temperature, 
size, distance and SED), the shell geometry (inner and outer 
radius, density variation), and the dust grains (species, optical 
depth at a reference wavelength, grain size distribution). Table Q] 
summarizes the model parameters that define the star+shell sys- 
tem. Several of these parameters are very sensitive to the mod- 
eling process and for these we have shown a best fit parameter, 
as well as a range of acceptable fits, in the table. Other param- 
eters are adopted based on reasonable assumptions because of 
limited choices, ancillary data or because the modeling process 
is not very sensitive to their choice. We describe the selection 
of these parameters and the estimation of the parameter range in 
the subsections below. 

3.1.1. Central Star 

We model th e central star using th e carbo n star model pho- 
tospheres of Ga utschv-Loidl et all (120041) . plac ing it at a 
distance of 50 kpc (LM C distance; see, e.g., iFeasl [19991 
Ivan Leeuwen et al" 120071) . We prefer these model spectra over 
a simple blackbody, as the inclusion of molecular absorption 
changes the photospheric features in the optical and NIR, which 
in turn affects the output spectrum (a 15-20% difference in IRAC 



fluxes between the 3000 K black body and the 3000 K mode l 
photosphere discussed below). The Gautschv-Lo idl etall ([2004) 
set includes 29 solar-mass, solar-metallicity models. As pho- 
tospheric molecular features are sensitive to metallicity (e.g., 
Mats uura et al.ll2002l) . we may not be able to faithfully repro- 
duce their strengt hs in LMC stars (Z/Z -0.3-0.5 for the LMC; 
IWesterlundll997h . but such differences for LPV 28579 would be 
overwhelmed by emission from its optically thick dust shell. The 
near-IR photometry suffers from circumstellar extinction, mak- 
ing it difficult to constrain the effective temperature and surface 
gravity. Therefore, we are only able to find a suitable model by 
using the source luminosity as a constraint, which we estimate 
to be (7690±400) L for the source using the Epoch 1 fluxes. 
Two models have luminosities of 6090 L , within 20% of our 
estimate for the star. They have effective stellar temperatures 
of 3000 K and 3100 K, and logg values of -0.55 and -0.49 
respectively. To model LPV 28579, we use the 3000 K photo- 
sphere. Keeping all other parameters fixed, we found that there 
is very little difference in the near- and mid-IR output SEDs if we 
picked the 3100 K model instead. The model photosphere has a 
C/O ratio of 1.3. While this value is typical for Galactic Miras, 
the C/O r atio may be higher for Magellanic Cloud carbon stars 
(see, e. g. . iMatsuura et alj|2005t IStanghellin i et alj|2005l) . which 
would change the optical depth due to increased molecular opac- 
ity. However, since we are limited by the small number of model 
photospheres available, we ignore these effects. 

3.1.2. Circumstellar Shell Geometry 

We make the simplifying assumption of spherical symmetry as 
we are more interested in constraining the dust properties in 
this study. We also assume a constant mass-loss rate and out- 
flow velocity y exp , leading to an inverse-square density distri- 
bution in the shell. The lower metallicity of the LMC suggests 
lower y exp values for the O-rich AGB outflows compared to their 
Galactic counterparts, but the out flow velocities of carbon stars 
are le ss sensitive to metallicity (lHabindll996l) . IWachter et alj 
(2008) compute carbon-star dynamical wind models for sub- 
solar metallicities. They find that the outflow velocities of solar 
metallicity models are about twice that of the LMC. These theo- 
retical predictions of lower outflow vel ocities at lower metal lic - 
ity are supported by the observations of lLagadec et alJ (l2010h . If 
the expansion velocity is assumed to be independent of distance 
in the shell, it has no effect on the shape of the output SED and 
cannot, therefore, be constrained by performing radiative trans- 
fer. For simplicity, we assume a value of w e xp=10 km s , which 
is in agreement with observations of Galactic C-ric h AGB stars 
of similar periods (see Fig. 4 of iMarigo et alJl2008l) . The mass- 
loss rate depends linearly on i/ exp (see Eq.[2]) for our simplifying 
set of assumptions, therefore the calculated rate will be lower if 

the outflow velocity is lower than the Gala ctic value. 

Radiative transfer codes like DUSTY dNenkova et al.ll2000h 
accept the dust temperature at the inner radius as input and cal- 
culate the inner radius on execution, and this temperature is typ- 
ically fixed at the condensation temperature of the dust species. 
This condition is not necessarily true for the highly evolved AGB 
stars as well as post-AGB stars. 2Du st adopts a more general 
approach (see lUeta & Meixnerll2003l) by using the inner radius 
to calculate the temperature structure in the shell. For sources 
with detached shells, the inner radius can be directly observed 
in high-resolution images. The shape of near- and mid-IR SED 
constrains the inner radius and therefore the temperature at this 
radius as well. The inner radius of the shell need not, therefore, 
correspond to the condensation radius, and the resulting temper- 
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ature may not be equal to the condensation temperature of the 
dust. 

The ratio R mdx /Ri n of the outer radius of the dust shell to 
the inner radius is kept fixed at 1000. The outer radius deter- 
mines the total shell mass, which in turn is related to the mass- 
loss timescale. While this timescale is an important quantity, we 
are more interested in obtaining the AGB mass-loss rate, which 
is only weakly sensitive to changes in the outer radius as long 
as it is large enough to include contributions to the flux at the 
longest wavelengths of interest (in our case, the 24 pm band ) 
from the coldest dust in the shell. Like Ivan Loon et al. (2003), 
we find that the output SED is insensitive to changes in R max 
out to ~30 pm. The mass-loss rate is, however, very sensitive to 
the value chosen for the inner radius. Observations of G alactic 
carbon stars (e.g., IRC+10216 TDanchi et all 19901 1 19951) a nd re- 
sults of radiative transfer modeling (e. g.. Ivan Loon et al.l[l999t 
iGroenewegen et ail 1998) suggest inner radii in the range ~2-20 
stellar radii (R*), while simple estimates based on energy bal- 
ance suggest that amorphous carbon dust should form within a 
few R„ (lHofnerll2007l) . We vary the inner radius within these 
bounds until the desired SED shape is obtained. 

3.1.3. Dust Grain Properties 

We model the dust around LPV 28579 using amor phous car- 
bon ( optical constants from the "ACAR" sample of Zubk o et al.l 
1 1996b and silicon carbide (optical constants from Pegourie 1988) 
grains in order to reproduce the strong 11.3 //m SiC feature and 
long-wavelength continuum observed. 

Ou r grain sizes are distributed according to the iKim et all 
(fl994h (KMH) prescriptiorfl, with a power-law falloff of index 
7=3.5. We vary the minimum grain size a m [ n and the scale height 
flo until the shape of the spectrum is reproduced. We find no- 
ticeable changes in the output SED only when a m j n is increased 
beyond ~0. 1 pm. This insensitivity to the small grain size limit 
is due in part to the fact that larger grains contribute more sig- 
nificantly to the dust mass, making it easier to place constraints 
on the scale height th an on the minimum grain size (a similar re- 
sult was obtained by ISpeck et al.ll2009l) . Given this distribution 
of sizes as input, 2Dust then internally calculates absorption and 
scattering cross-sections and asymmetry factors for an "average" 
spherical dust grain from the Mie theory dMid fl908) assuming 
isotr opic scattering. We run the code in the Harrington averag- 
ing (lHarrington et al.l[l988l) mode, which gives an average grain 
size of ~0.1 pm, the typical value used in sin gle-size models 
(e.g.. Ivan Loon et aT]|1999l ; lGroenewegenl2006l) . 

2Dust requires as input the optical depth at a reference wave- 
length, which we choose to lie near the center of the SiC emis- 
sion (1 1.3 pm). We iteratively adjust this optical depth and rela- 
tive abundance of SiC until the shape of the feature is satisfacto- 
rily reproduced. 

3.1.4. Fitting procedure 

2Dust uses the dust opacities, dust grain and shell properties and 
the density variation as input to first calculate the density pj n at 
the inner radius from the relation 

ta = Ki J p(r)dr (1) 

where is the optical depth at the reference wavelength A, and 
Kx is the opacity averaged over dust size and composition. The 

7 n(a) oc aT y exp [-a/ao] with a > a m i n 



total dust mass in the shell is calculated once p m is known. This 
is divided by the wind crossing time in the shell to obtain 

M = 4n—R in v exp (2) 
ka 

The dust mass-loss rate can be calculated once a value for the 
expansion velocity is chosen. 

The overall fitting procedure is summarized as follows. 
Starting with the model photosphere that comes closest in flux 
to the luminosity of the source (see Sect. 13.1.11 ) we modify the 
inner radius, optical depth, SiC abundance and the range of grain 
sizes until the observed SED as well as the SiC feature strength is 
reproduced. In order to compare the results of 2Dust with our ob- 
servations, we generate synthetic photometry in the UBVIJHK^ 
as well as IRAC and MIPS24 bands in a manner similar to 
ISrinivasan et al.l (120091) . At each step, we scale the 2Dust out- 
put SED to the Epoch 1 8 pm flux for comparison. The scaling 
affects the luminosity as well as the mass-loss rate, and in Sect. 
13.1.51 we investigate the quality of the fit by varying this scale, 
which results in a range of predicted luminosities and mass-loss 
rates. 

Starting with the r e ff=3000 K photosphere, we first increase 
the optical depth until we roughly match the optical extinction as 
well as the mid-IR continuum observed in the photometry. The 
inner radius is now varied to improve the fit in the near-IR and 
IRAC part of the SED (the outer radius is kept fixed at 1000 R in , 
see Sect. l3.1.2l >. At this point, we compare the SiC features in the 
model and the IRS spectrum, and change the relative abundance 
of SiC in the model if required. In general, this change affects 
the IRAC fluxes, and we compensate for this by modifying the 
optical depth and/or inner radius as needed. Finally, we vary the 
minimum and maximum grain sizes to improve the fit. As noted 
in Sect. 13.1.31 changes in the minimum grain size do not signif- 
icantly affect the spectrum. We terminate the iterative procedure 
described above when we are able to satisfactorily match the op- 
tical, near-IR and IRAC/MIPS24 photometry for LPV 28579, as 
well as reproduce the 11.3 pm feature in its IRS spectrum. The 
resulting best-fit model is described below. 

3.1.5. 2Dust Fit Results and Sensitivity To Parameter 
Variation 

Table [TJ summarizes the fixed parameters as well as the values 
of the free parameters corresponding to the best-fit model. The 
2Dust output spectrum is superimposed on the photometry and 
IRS spectrum in Fig. [3] The luminosity of the model photosphere 
is lesser than the value calculated for LPV 28579 from photom- 
etry, so the entire SED has to be scaled up. The value of this 
scale factor decides the predicted luminosity and mass-loss rate, 
as well as the fit quality. We find that we are able to fit the pho- 
tometry as well as the IRS spectrum for a scale factor of 1.08, 
which corresponds to a luminosity of about 6580 L Q . Taking into 
account the photometric errors, a 1% to 15% raise in model lu- 
minosity (corresponding to luminosities in the range 6150-7010 
L ) result in good fits (;$ 7% variation around the best-fit value). 
To demonstrate that this uncertainty is small compared to the 
difference in luminosities between the two epochs, we scale the 
best-fit model down to fit the Epoch 2 photometry and find a 
luminosity of about 4810 L Q (~27% fainter). The scaled best- 
fit dust mass-loss rate for Epoch 1 data is found to be about 
2.5 x 10~ 9 M yr" 1 . The optical depth at 11.3 pm is 0.28. The 
inner radius of the dust shell is 4.4 R, (1430 R ) and the temper- 
ature at this distance from the star is 1310 K. 
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The results of the model fitting are most sensitive to changes 
in R m , the optical depth, the SiC dust abundance and the value 
of oo- When we have created a grid of models, we will be able 
to derive chi-squared errors for the parameters. For the moment, 
we determine the acceptable range of variation in the parame- 
ters mentioned above as follows. We create boundary ranges on 
the data by allowing the photometry to change by up to 3 times 
the uncertainty and by restricting the spectroscopy to a range 
bounded by its lcr uncertainty, because of lower signal-to-noise 
ratio associated with the spectra. Using this range of data vari- 
ation, we determine by eye the range of variation in the model 
fit due to changes in a certain parameter while keeping all other 
parameter values fixed at their best-fit model values. The range 
of parameters thus calculated are shown in TableQ] Fig.|4]illus- 
trates this process by showing the range in acceptable model fits 
as a result of changing the optical depth. 

The process described above helps us place a rough con- 
straint on the variation of the important output parameters - the 
luminosity, mass-loss rate and the dust temperature on the in- 
ner radius. We find that changing the inner radius of the dust 
shell causes the largest variation in both the temperature at the 
inner radius as well as the mass-loss rate. We find T m values be- 
tween 1260 K and 1380 K. The mass-loss rate is in the range 
(2.4 - 2.9) X 10 -9 M yr _1 . The variation of parameter values 
thus gives us an estimate of the required data quality if the de- 
rived mass-loss rates are to be within ~15%, assuming our value 
for fexn is correct. 



3.2. Modeling The Circumstellar Gas 

The moderate optical depth of the circumstellar shell should 
be effective in attenuating photospheric features. However, LPV 
28579 has strong absorption features around 13.7 fim. The fea- 
tures therefore originate in the circumste llar region (see, e.g., 
Ziilst ra et alJ 120061: iMatsuura et alJ 120061) . In Galactic carbon 
stars, there is also a contribution to this wavelength range due 
to HCN, but due to the low ab undance of nitrogen, this effect is 
not pronounced in LMC stars (Matsuura et alJl2006l) . 

We simulate the molecular ban ds of CO and C 2 H2 in LPV 
28579 using the slab model code of Matsu ura et alJ (|2002) with 
recent line lists published bv lRothman et alJ (120091) as part of the 
H1TRAN database. Assuming that the majority of the molecu- 
lar absorption comes from the circumstellar shell, the excitation 
temperatures are chosen to be cooler than the temperature at the 
inner shell radius, i.e., T ex < T[ n . Based on these considerations, 
we derive excitation temperatures of 1000 K for both species, 
with column densities of 10 19 cirT 2 and 3 x 10 21 crrT 2 for C2H2 
and CO respectively. Fig. [2] shows the locations of the spectral 
features for the parameters above. The resulting (gas+dust) fit 
is in good agreement with the spectrum and is plotted in Fig. 
[3] The CO band is at the edge of the IRS spectrum so we do 
not have observational constraints on it, but we are able to re- 
produce the 13.7 fim feature using C2H2. We note here that the 
model is useful only for identifying the molecular bands and 
to roughly estimate the excitation temperatures and densities of 
these molecular species, it can not be used to determine the de- 
tails of the atmospheric structure. A determination of the mass- 
loss rate would also require knowledge of the location of the 
slab of C2H2 in the circumstellar envelope. As this quantity is 
very uncertain, we do not calculate a gas mass-loss rate using 
this method. 



4. Discussion 

The long-term goal of this study is to use 2Dust modeling of car- 
bon star circumstellar shells to fit the photometry of the ~7 000 
C-rich SAGE AGB candidates, along with any carbon stars 
found among the ~ 1 400 sources in the SAGE extreme AGB list, 
in order to derive their mass-loss rates and to estimate the total 
mass-loss return from AGB stars to the LMC ISM. We are also 
deriving the dust properties for O-rich AGB stars in a parallel 
study (Sargent et al. 2010, in press). In this section, we weight 
the simplifying assumptions used in our modeling of LPV 28579 
against this future goal. Our study of the variation in model pa- 
rameters demonstrates that we place reasonable constraints on 
the inner shell radius, the temperature at this radius, dust grain 
properties and the resulting mass-loss rates. We discuss the re- 
sults of our modeling in detail below. 

4.1. Dust Properties 

Figure [3] illustrates the good agreement between our model us- 
ing an AmC+SiC dust composition and the overall continuum 
as well as the 11.3 fim feature observed in the spectrum. Our 
combined gas+dust model improves on the 2Dust prediction by 
fitting the 13.7 /mi absorption. We note two minor discrepancies: 
the silicon carbide feature seems broader than the model predic- 
tion, and the model fit systematically overestimates the flux at 
wavelengths longward of ~20 /urn. A broade ning of the SiC fea- 
ture may indicate self-absorption (see, e.g.. ISpeck et al.l["l997h 
corresponding to a high optical depth. The long-wavelength dis- 
agreement is probably due to the poor background subtraction 
(see Sect. \2.3\ rather than due to a modeling issue. Given the 
low signal-to-noise of the spectrum, it is not certain if these dis- 
crepancies are real; we are therefore unable to justify a detailed 
refinement of our RT model. 

Our best-fit model predicts a temperature of 1310 K at the 
inner radius of the dust shell, which is located at 4.4 R„. In Sect. 
13.1.51 we estimate the inner edge to lie within (4.0^4.8) R t (cor- 
responding to temperatures of 1260-13 80 K), which is consis- 
tent with the theoretical predictions of iHofnerl (120071) for the 
condensation of amorphous carbon dust (7\. om /=1500 K, ^; n =3 
R*). The range of temperatures are slightly warmer than would 
be expected for LMC carbon stars of comparable mass-loss rate s 
(see, e.g., Fig. 4 and Equation 1 of iGroenewegen et alj |2009). 
This may be due to the fact that we are using a warmer central 
star. Here we are limited by the availability of the colder pho- 
tospheres of comparable luminosity. Moreover, as already dis- 
cussed, the moderate optical depth prevents us from using the 
observed near-IR colors to estimate the central star's effective 
temperature. 

The KMH grain size distribution for the best-fit model is de- 
fined by y=3.5, a m i n =0.01 //m and ao=1.0 fim. As noted in Sect. 
13.1.31 the shape of the output SED is sensitive to the size of the 
largest grains, but does not vary appreciably with changes in the 
minimum grain size up ~0.1 /im. About 15-18% of the mass is 
contained in grains with diameters exceeding 1 /zm, depending 
on the minimum grain size. For comparison, only 0.1% of the 
dust mass ejected from IRC+10216 (L ~1.5 x 10 4 L , i/ exp =15 
km s ~', dust M LR ~8 x 10 M Q yr~') is in micron-sized parti- 
cles (lJuralll994l) . 

The SiC content derived from the best-fit model is com- 
parable to values derive d for dusty Galactic carbon stars 
(IGroenewegen et aTlll998l) . The lower metallicity of the LMC 
would imply a lower s ilicon abundance in general, and 
IGroenewegen et all d2007l) find lower SiC/AmC ratios for the 
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LMC compared to Galactic stars. By studying the effect of vary- 
ing the SiC abundance on the output SED, we note that the SiC 
content can be in the range 10-16%. In this paper, we assume 
a gas:dust ratio ¥=200. This valu e was dete rmined for carbon 
stars in the solar neighbourhood dJur a 1986). An accurate de- 
termination of the gas:dust ratio is still to be made for LMC 
carbon stars. Some observational and the oretical studies (e.g., 
iMatsuura et alj|2005t IWachter et"ai]|2008l) indicate that the car- 
bon star gas:du st ratio may not be metallic i ty dependent. In fact, 
recent studies dGroeneweeen et al.l l2007t [Wachter et al. 2008; 
iMattsson et al.ll2008l) have found that the mass-loss rates of car- 
bon stars may not be metallicity dependent. The gas:dust ratio 
can be estimated from the mass fraction of SiC in the dust if we 
also know the fraction of the total Si mass that is in the dust. 
Assuming that the Si abundance scales with metallicity (The 
LMC silicon abundance is n ot well constrained; see the discus- 
sion in Matsuura et al. 2005) we obtain ¥ in the range 200-570. 
However, we have assumed here that all of the Si goes into the 
dust. In reality, the fraction of Si in dust depends indirectly on 
the amount of carbon in the dust. For example, if half of the car- 
bon remaining after CO formation goes into dust, then 8-13% 
of the Si is in dusQ The situation improves for higher C/O ratio 
- for C/0=1.8, the Si dust fraction is 20-34 %. An alternative 
way to estimate the gas: dust ratio is to use the Ivan Loonl (120001) 
relationship ¥ ~ Z~ 1±0 - 3 for carbon stars, which corresponds to 
¥=300-1000 for the LMC. Due to the lack of constraints on 
the various quantities that determine the gas:dust ratio, we use 
¥=200 when quoting a total mass-loss rate for LPV 28579, not- 
ing that the mass-loss rate could be up to about 5 times larger 
based on the calculations above. 



4.2. Mass-loss Rate 

Our 2Dust model predicts a dust mass-loss rate of 2.5 x 10~ 9 
M yr -1 . For a gas:dust ratio of 200, the total mass-loss rate 
is then 5 x 10~ 7 Mp yr" . This rate is on the low end of the 
Ivan Loon et al. I ( fl999h mass-loss rates for the brightest, most 
obs cured LMC carbon stars a s well as the rates calculated for 
the iGroenewegen et al.1 d2009) sample. The mass-loss rates cal- 
culated from the [3.6]-[8.0] and K-\S.O] co l ors of LPV 28579 
using equations 1 and 2 of IMatsuura et al.l (120091) are slightly 
higher (1.4 x 10~ 6 M yr" 1 and 2.5 x 10" 6 M yr" 1 respec- 
tively). However, our dust MLR is consistent with the value de- 
termined from the 8 fim flux - using the excess-MLR relation 
for extreme AGB stars from Paper I, we get a rate of 3.9 x 10~ y 
M yr . Our total rate agrees with the values calculate d us- 
ing the period-MLR relatio ns of|Vassiliadis & Wood! (1 19931) and 
IGroenewegen et al l (1 19981) (10~ 7 M yr" 1 and 7.4 x 10" 7 M 
yr~' respectively, for a period of 356. 2d). The value obtained for 
the mass-loss rate is sensitive to the v exp adopted (see Eq. [2]). 
The outflow velocities for Galactic carbon stars a re typically 10 
km s~ 1 , but can be up to about twice as high. The lWachter et alj 
(2008) models predict LMC outflow velocities to be about half 
that of comparable Galactic sources. This range of velocities in- 
troduces an uncertainty of a factor of two in the calculated mass- 
loss rate. 



8 Here, we have assumed that carbon and silicon are the only con- 
stituents of the dust. A more realistic estimate would include metallic 
iron dust, which somewhat counteracts the low Si fraction and lowers 
the gas:dust ratio. 



4.3. Evolutionary State 

In order to understand our L and r e jf values for LPV 28579 
in the context of its evolutionary state, we identify its loca- 
tion hi_Jhe_Iiertzsgrung-Russel (HR) diagram in relation to 
the iMarigo et al.l (|2008) theoretical isochrones. We choose a 
metallicity of Z=0.006, which is consistent with observational 
estimates for the metallicity of the LMC. We generate these 
isochrones using their CMD interactive web interface^ so that 
we can also extract the initial stellar masses for the stellar tracks. 
Figure [5] shows the location of LPV 28579 on the HR diagram. 
Also shown are selected isochrones for Z=0.006 (ages ranging 
from log t /Gyr=9.0 to 9.8). We find that the log t /Gyr=9. 1 to 9.6 
isochrones are able to produce the range of luminosities we have 
estimated for LPV 28579 in this stud y, which gives us a n age 
in the range 1.6-4 Gyr. As noted bv IMarigo et alJ (|2008), this 
is the typical age of the O-rich AGB population that undergoes 
a smooth transition to C-rich chemistry due to third dredge-up 
(also see the discussion of the multiple O-rich populations in 
Paper I). The initial masses correspo nding to these isochrones 
were i n the range 1 .3-2 M . Figure 3 of lGroenewegen & de Jond 
(11994 shows the evolution of 1.25 and 5 M AGB stars in 
period-luminosity space. Based on its luminosity and primary 
period, LPV 28579 falls slightly upward of the 1.25 M track, 
providing a mass comparable with our estimate using theoretical 
isochrones. 

The derived age (and hence progenitor mass) range is depen- 
dent on the metallicity of the isochrones as well as our T e ff esti- 
mate for LPV 28579. Figure[5]demonstrates the effect of chang- 
ing the metallicity (compare the log f/Gyr=0.94 isochrones for 
Z=0.006 and 0.008). At lower metallicities, the isochrones pre- 
dict carbon stars at lower luminosities. The optically thick shell 
of LPV 28579 makes the precise determination of T e s difficult 
(see 13.1.1b so that our choice of model photosphere is made 
on the basis of the luminosity estimated from photometry. The 
range of ages and masses from isochrone comparison is depen- 
dent on how well we can constrain the effective temperature. If 
we assume that T e ff is cooler than 3100 K, we obtain the age and 
mass range above. For r e ff under 3050 K, the range reduces to 
log f/Gyr=9.2-9.5 giving masses in the range 1.4-1.7 M . 

The total mass-loss rate of 5.0 x 10~ 7 M yr~' is comparable 
to the nuclear-burning rat e at the luminosity of the star (see Fig. 
9 in lvan Loon et alj[l999l) . which suggests that LPV 28579 may 
undergo the superwind phase at the end of its life on the AGB. 

5. Conclusions 

We have modeled the circumstellar shell of the carbon-rich 
asymptotic giant branch star LPV 28579 using the radiative 
transfer code 2Dust to fit the SAGE photometry and S AGE-Spec 
spectroscopy available for the source. Using a mixture of amor- 
phous carbon and silicon carbide dust with 12% SiC by mass, 
and dust grains ranging in size from 0.01 /mi to ~ 1.0 fim, we find 
a luminosity of 6580 L and a dust mass-loss rate of 2.5 x 10~ 9 
M yr -1 . The inner radius of the shell is at 4.4 times the stel- 
lar radius and at a temperature of 1310 K. The optical depth at 
11.3 fim is found to be 0.28. We model the molecular features 
in the circumstellar shell with slab models of C2H2 and CO, and 
calculate column densities of 10 19 cirr 2 and 3 x 10 21 cm 2 for 
the two species respectively. The model SED when combined 
with the gas model shows excellent agreement with the observed 
spectrum and photometry. The total mass-loss rate assuming a 

9 http://stev.oapd.inaf.it/cmd 
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gas:dust ratio of 200 is comparable to estimates derived from 
empirical relations involving the period, color and mid-infrared 
flux. LPV 28579 has properties typical of dust-enshrouded LMC 
carbon Miras, and the model derived in this work defines a base- 
line set of carbon star dust properties in our upcoming model 
grid. 
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Table 1. 2Dust Parameters For LPV 28579 



niuiuapncrc muuci 




Teff(K) 


1 AAA 

3000 


i*(L ) 


iff OA 1 C A 1A1 A\/? /T" 1 

6580 (6150-7010) (Ep 




4olO (tpocn Z) 




jzj (tspocn i ) 


c/o 


1.3 


Dust shell properties 






4.4 (4.0-4.8) 


^out(^in) 


1000 


Density profile 


p(r) ~ r~ 2 


VexpCkmS" 1 ) 


10 


Gas:dust ratio 


200 


Dust grain properties 


Species 


AmC e +SiC rf 


SiC fraction 


12% (10%-16%) 


t(11 yum) 


0.27 (0.25 - 0.275) 


Size distribution 


KMH e 




a^Qim) = 0.01 (<0.1) 




flo(A^m)= 1 (0.75-1.3) 




y = 3.5 



Mass-loss rate and dust temperature 

M^xlO- 9 M yr 1 ) 2.5 (2.4-2.9) 
M^xlO" 7 M yr 1 ) 5.0(4.8-5.8) 
T in (K) 1310(1260- 1380) 



Notes. (a) Photosphere model from Gautschy-Loidl et al. (2004j). m Un certainty in parameter value (see Sect. l3.1.31 for detail s). M Amorphou s 
carbon grains, p = 1.8 g cm -3 , optical con stants fromlZubko et all d 19961) . Id> SiC grains, p = 3.22 g cm 4 , optical constants from Pegourie (1988). 
w Size distribution from lKim et all j!994l) : 

n(a) ~ a~ y exp (—a/ag) with a > a min . 
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\ (/jm) 

Fig. 1. A compilation of photometric and spectroscopic infor- 
mation for LPV 28579. The solid circle and error bar show the 
mean OGLE / band flux and its range of variation respectively. 
The 2MASS (open triangles), IRSF (filled triangles) and SAGE 
(Epoch 1 : open squares, epoch 2: filled squares) fluxes are also 
shown. The solid curve is the SAGE-Spec IRS spectrum. 
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Fig. 2. The synthetic spectra for CO (dashed line, column density 
= 3xl0 21 cm 2 ) andC2H2 (thin solid line , column density = 10| 9 
cm" 2 ) from the HITRAN database and iRothman et alJ d2009). 
The combined spectrum (thick solid line) used to model the 
molecular contribution to the observed spectrum is also shown. 
The spectra for CO and C2H2 are scaled down for clarity. 



11 



S. Srinivasan et al.: Modeling LMC carbon stars 




Fig. 3. The results of modeling the dust and gas around LPV 
28579. The SAGE Epoch 1 (red dots) and epoch 2 (black 
dots) photometry as well as SAGE-Spec spectrum (black curve) 
shown with the best-fit 2Dust RT model (green curves) scaled to 
fit Epoch 1 and Epoch 2 photometry . The photosphere model 
corresponding to the best-fit SED is also shown (dashed curve). 
The model SED is convolved with the synthetic molecular spec- 
trum shown in Fig. [2] to reproduce the 13.7 jjm feature (pink 
curve). 
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Fig. 4. The sensitivity of the 2Dust model output to optical depth. 
The SAGE Epoch 1 photometry forLPV 28579 is allowed a 3cr 
variation (magenta error bars), and the spectroscopy is kept to 
within lcr (gray curves). Model SEDs for 11.3 fim optical depth 
values spanning this range of allowed fluxes are also shown. The 
optical depth for good fits ranges from 0.26 (blue curve) to 0.29 
(red curve). The best-fit curve for r=0.28 (green) is also shown 
for comparison. 
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Fig. 5. LPV 28579 (black dot with error bar representing range 
of calculated Epoch 1 luminosities) on the L - T e ff diagram. 
The solid and dot-dash ed lines, tagge d according to their age 
(log f/Gyr), are the lMarigo et alj 02008) isochrones forZ=0.006. 
For comparison, the Z=0.008 track for logf/Gyr=9.4 is also 
shown (dashed line). 
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